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The non-equilibrium distribution of colloids in a polymer solution under a temperature gradi- 
ent is studied experimentally. A slight increase of local temperature by a focused laser drives the 
colloids towards the hot region, resulting in the trapping of the colloids irrespective of their own 
thermophoretic properties. An amplification of the trapped colloid density with the polymer con- 
centration is measured, and is quantitatively explained by hydrodynamic theory. The origin of the 
attraction is a migration of colloids driven by a non-uniform polymer distribution sustained by the 
polymer's thermophoresis. These results show how to control thermophoretic properties of colloids. 

PACS numbers: 82.70.Dd, 66.10.Cb, 47.57.J- 



Introduction. — Gradients of thermodynamic variables 
such as temperature, chemical potential, and osmotic 
pressure cause migration of molecules and small particles 
in fluids pQ. For example in biological cell, coupling of two 
gradients is often used to promote molecular transport 
against one of the gradients as in chemiosmosis[2j. In 
physics and chemistry, novel methods to utilize phoretic 
properties (electro-, thermo-, and diffusiophoresis) are 
proposed for transporting and screening particles in lab- 
on-chip [3j [4] or designing active materials [5 j. In ther- 
mophoresis, the speed and direction of migration along 
a temperature gradient are characterized by Soret coef- 
ficient, which is generally material-dependent [6j 13 [HJ |9j 
[TO] . Diffusiophoresis is a similar phenomenon where col- 
loids move along a gradient of solute molecules [HI [12] 
[13]. The benefits of them have been used for efficient 
and amplified transport in microfuidics[T4 ] [T5 ] [T6] . How- 
ever, strong material dependence of Soret coefficients, 
and difficulty in keeping a stationary solute gradient of 
diffusiophoresis prevent further development of applica- 
tion of nonequilibrium force by the phoretic effects. Un- 
derstanding physical mechanism of thermophoresis has 
made marked progress recently [T7| fT8| [19], however, it 
is still challenging to control magnitude of the force to 
overcome these problems. 

In this paper, we demonstrate that a suitable coupling 
of thermophoresis for polymer solute molecules and dif- 
fusiophoresis for a target particle (colloid) largely resolve 
the problems. More specifically, we report experimental 
and theoretical studies on a phoretic motion of colloidal 
particles in a polymer solution under a temperature gra- 
dient. We find that a Soret coefficient of a colloid is 
sensitive to the polymer (PEG) concentration; as increas- 
ing the amount of polymer, the Soret coefficient reverses 
its sign and its magnitude outweighs by far its intrin- 
sic value at the highest polymer concentration studied. 
The dependence of the Soret coefficient on the polymer 
is experimentally determined and is corrobolated by our 
hydrodynamic calculations. This novel effect allows us 
to transport and trap colloids at any desired position by 
suitably controlling a temperature distribution and the 
polymer concentration. 



Experiment. — In a thin chamber containing a solution, 
a steep temperature gradient up to 1 [K//im] was cre- 
ated while keeping the local temperature rise quite small 
(AT - 4K) by focusing an infrared laser (NdiYAG, 1064 
nm, power ~4mW before an objective lens) on a light 
absorbing metal coated surface of the thin glass cham- 
ber (Setup I) as shown in Fig. la [20]. Using this setup, 
we measured the distributions of polystyrene beads. Fig- 
ure lb shows fluorescence intensity of 100 nm diameter 
beads (F8803; Invitrogen) in a water solution. The fluo- 
rescence intensity became lower at the hot region around 
laser focus, indicating polystyrene beads escaped from 
the hot region due to thermophoresis. The effects of laser 
trapping and convection are negligible compared with 
thermophoresis of the beads in this setup. Polystyrene 
beads and typical biomolecules such as DNA migrate to 
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FIG. 1: a) (color online) Setup I: The laser fed into a thin 
glass chamber (2/xm) is focused on the top surface coated by 
a thin Cr layer for light absorption, b) Distribution of 100 
nm fluorescent polystyrene beads in water, c) Distribution 
of 100 nm fluorescent polystyrene beads in PEG solution, d) 
Setup II: The laser fed into a chamber (20/xm thick) is used 
for directly aqueous solution heating, e) Cluster of 500nm 
polystyrene beads in PEG solution. Scale bar: 10/xm 
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FIG. 2: (color online) a) Profile of temperature (green: rect- 
angle, right ordinate) and density of fluorescent PEG (blue: 
circle, left ordinate) as a function of radius from the center 
of laser focus, b) Density of fluorescent PEG as a function of 
temperature increment. 

colder regions at room temperature 9 . However, when 
we added a small amount of neutral polymer, polyethy- 
lene glycol (PEG6000, MW7500, 3.5wt%) in lOmM Tris 
buffer, beads and DNA molecules [21] migrated and were 
trapped at the hot region regardless of the sign of ther- 
mophoresis (Fig.lc). The range of attraction reached 5 
to 10/im, corresponding to the size of the temperature 
distribution (see Fig. 2a), exceeding by far about 1/im in 
optical tweezers. As we moved the laser or chamber at 
several micrometers per second, the trapped polystyrene 
particles moved with the hot region (see the support- 
ing information [22]). This indicated that the sign and 
the magnitude of thermophoresis were modified in the 
presence of polymer under a temperature gradient. This 
effect was also observed in a bulk heating configura- 
tion as shown in Fig. Id (Setup II). Direct absorption of 
laser light with longer wavelength(1480nm, 25mW, Fu- 
rukawa) in the water produced a local temperature gradi- 
ent, which creates a 3D colloidal aggregation within the 
hot region (Fig. 2e) 

Polymer Distribution. — To elucidate the mechanism 
of the attraction, we measured a profile of temperature 
and density distribution of PEG near the laser focus in 
the Setup I. The temperature profile and the density of 
PEG molecules in chamber were shown in Fig.2a[23]. In 
this setup, temperature profile and density profile were 
axisymmetric. Thus the data shown below are averaged 
along azimuthal direction. The temperature in the laser 
heating spot increased approximately 1 K with a 1 mW 
increase in the laser power. We found that the PEG 
molecules were depleted by about 15% due to PEG ther- 
mophoresis when the temperature increment was about 
3K in the center (Fig. 2a, b). The depletion of PEG can 
be described as a steady state distribution by the bal- 
ance between thermodiffusion and diffusion as follows. 
The flux of PEG molecules obeys the relation, 

j p = -D p Vc - cD p T VT, (1) 

where c 1 j pi D p 1 D P ^ are the polymer density, flux, dif- 
fusion constant, and thermodiffusion constant, respec- 
tively. In the steady state (j p = 0), the density satisfies 

dc/c = -{D p T /D p )dT = -S p T dT, (2) 



where the Soret coefficient defined by Sj, = D p T jD p is 
measured by S V T = -(l/c)(dc/dT) ~ -A In c/AT . We 
obtain S p = 0.046 ± 0.005 [K' 1 } for PEG5000 from the 
slope in Fig. 2(b), in agreement with ref . [24] . 

Trapping Ability. — Next we examined how the accu- 
mulation of beads is dependent on the concentration of 
polymer and on the local temperature increase by mea- 
suring the density distribution of beads. The fluorescence 
images of beads were recorded for different PEG6000 con- 
centrations. In Fig. 3, spatial distributions of fluorescent 
polystyrene beads (100 nm, 0.05%) were displayed. Each 
image was integrated for 300 frames in 20 seconds and 
was linearly scaled to 256 grey levels. In 0% solution, 
the density decreased from the outer region to the center 
(Fig. 3a). This depletion of beads is explained quanti- 
tatively by thermophoresis, as described below. In 1% 
solution (Fig. 3b), there was still no trapping, but the 
degree of depletion is diminished compared to the 0% so- 
lution. In 2% solution, the trapping phenomenon began 
to appear and the fluorescence intensity roughly doubled 
in the hot region (Fig. 3c). In the 3.5% and 5% solu- 
tions, the density increased approximately 10- and 100- 
fold at the center compared with the surrounding con- 
centrations, respectively (Fig. 3d, e). 

In Fig. 4, the particle density at each distance from the 
center is plotted against temperature at the same radius 
(Fig. 4). The result for 0% solution can be described as 
a steady state distribution due to thermophoresis. The 
density of beads, n, obeys Eq. ([2| with replacement of c 
by n and S J by 5^, where is the Soret coefficient of 
the beads. Neglecting temperature dependence of 5^, an 
exponential distribution is obtained [14 

n = n exp[-4(T-T )]. (3) 

The slope of the curve in Fig. 4 gives the Soret coefficient 
as Sj, = 0.35 ± 0.03 [K _1 ] in agreement with previous 
results [9 . Importantly, even in the presence of polymer, 
the particle density varied exponentially as a function of 
the temperature increment (see Fig. 4). Therefore a slope 
of each curves gives an effective Soret coefficient, S^. We 
plot S?p as a function of PEG concentration in the inset 
of Fig. 4. As a first-order approximation, decreased 
linearly with the increase of polymer concentration ac- 
cording to the relation 

S* T = S b T - X c (4) 

with x — 35.4 ± 2.0 [K _1 ] as the best fit, and where c 
is the concentration of PEG in weight fraction, or x = 
441±25[nm 3 K _1 ] for c in volume density. Hence, the sign 
and the magnitude of the Soret coefficient was controlled 
by varying the concentration of polymer. 

Theory. — The main results obtained in our experi- 
ment are quantitatively explained by the hydrodynamic 
theories [T2] [T8] [T9] . We consider a spherical particle of 
radius a in a polymer solution, see the inset of Fig. [5] In 
a dilute regime valid to our experiments, PEG distribu- 
tion around the particle obeys Boltzmann distribution, 
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FIG. 3: Fluorescence intensity of 100 nm beads under temperature gradient for different polymer concentrations, a) 0%, b) 
1%, c) 2%, d) 3.5%, and e) 5% PEG 6000 solutions, respectively. The laser is focused at the center. The scale bar: 10/xm. 
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FIG. 4: (color online) The relationship between the temper- 
ature increment and the density of the beads divided by its 
background value, for different PEG concentrations indicated. 
Filled symbols are the experimental data obtained from Fig. 3. 
The solid lines are the fits with Eq. ([3]), from which the effec- 
tive Soret coefficients, S^, were obtained. Inset: is plotted 
as a function of the PEG concentration in weight fraction, c. 
The solid line is the fitting with Eq. (7]). 



c = Co exp(— U{r) jk B T), where Co is the concentration 
at infinity, k B is the Boltzmann constant, and U(r) is a 
short range potential for PEG molecules. As a PEG is 
non- ionic and inactive, U represents an entropic repul- 
sion from the colloid surface with its interaction distance 
A. Note that A also defines a length scale of the deple- 
tion interaction between colloids at equilibrium [25} [26]. 
In a temperature gradient dT/dz applied along z direc- 
tion, a PEG gradient dc^/dz is developed according to 
Eq. As Co changes only gradually at a scale of A, the 
steady state distribution is given by c ^ co(z)e~ u ^ r ^ kBT . 
Neglecting interactions between polymers, an osmotic 
force density is given by / = — c(r)VC/(r), which is 
written in, / = V{k B T5c) + k B T(S? - ^)£cVT, where 
Sc = c-c = c (e- c/ / /esT -l), and Vc = -c^S^VT was 
used. The velocity field around the colloidal sphere, v, is 
obtained by solving the Stokes equation, rf\7 2 v = Vp — /, 
with the incompressibility condition, V-u = 0, where n is 
the viscosity of the solution and p is the hydrostatic pres- 
sure. The Stokes equation is rewritten in the following 
transparent form, 



where / (r) = k B Tc(S^_- 1/T)(1 - e~ u ^ BT )dT/dz and 
c = co(z = 0). In Eq. (J5|, the osmotic pressure, k B T5c, is 
compensated by a hydrostatic pressure p. (This ensures 
the absence of osmotic force on the particle proportional 
to its volume, ~ k B Ta 3 5c.) The nonequilibrium force, 
fc,e z in Eq. (J5|, on the other hand, is balanced with the 
shear viscous force, rjS7 2 v, leading to a fluid flow relative 
to the colloid surface and a phoretic motion of the colloid 
at a velocity u relative to the fluid at infinity. Solving 
Eq. ([5| in the colloid-reference frame with the no-slip 
boundary condition, and summing up the resulting fluid 
stress over the colloidal surface, we obtain the total force 
acting on the colloid as F = —Qnnau + 47re 2 J a °° r(r — 
a)/nrfr[I5]. This must be zero because no external force 
is acting on the colloid. Choosing a specific potential of 
hard-core type (i.e., U = oofora<r<a + A and zero 
otherwise), we finally arrive at 



k B T 
3 V 



(S P T 



)A 2 cVT. 



(6) 



Plugging then Eq. §6y into an expression of the density 
current of colloids, J = nu — DVn — uDtVT[19 , and 
rewriting it in a form J = —DVn — nDS^VT, we see 
that an effective Soret coefficient is given by 



UT 1 UT 1 



2tt(5P - ha\ 2 c 



(7) 



v V 2 v = V(p - k B TSc) + f (r)e z 



(5) 



which confirms the linear dependence on c in Eq. @, 
where x = 2tt(S£ - l/T)aA 2 . 

Equation (JTl) is compared to our experimental data 
in Fig. [4] inset by employing = 0.056 for PEG with 
MW=7500 estimeted from [24 . The best fit gives A ~ 5.2 
nm, which is close to the PEG gyration radius ~ 3 
nm [27]. We consider the agreement satisfactory, as A 
should be of the order of a PEG size. A further consis- 
tency check was made by looking at the dependence of 
S?r on the particle radius a. To focus only on this depen- 
dence of we define an effective depletion volume at 
nonequilibrium V = x/(Sf. — 1/T). First, we performed 
the series of experiments for four different colloid sizes, 
and extracted V from the data on S^, which are plotted 
in Fig.[4]as a function of the diameter 2a. The data shows 
an overall linear dependence, and the agreement with the 
theoretical prediction for no-slip condition, V = 27rA 2 a, 
obtained from Eq. ^ is good. Note that we used A = 5.2 
nm obtained before, and thus no adjustable parameter is 
assumed. Nevertheless, it would be important to examine 
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FIG. 5: (color online) Particle diameter dependence of the 
depletion volume V. Dotted line and solid line represent the- 
oretical predictions for non-slip and Navier boundary condi- 
tion, respectively. Inset : (top) schematic illustration of a 
colloid particle in polymers solutuion. (bottom) Soret coeffi- 
cients of beads for different sizes. 

effects of slip of the polymer solution upon the enhance- 
ment of the Soret effect, since the perfect slip condition 
predicts a different scaling, V ~ Aa 2 [19]. A bound- 
ary condition for a general slip case is a so-called Navier 
condition [28 j, given by a t = rjv t /£ at r = a, where a t 
and vt are the tangential components of the fluid stress 



and velocity. The parameter £ defines a slip length on 
the surface (the no-slip or perfect slip limit is attained 
for £ = or £ = oo)[29 . In this Navier case we ob- 
tain V ~ 27raA[aA + 2£(a + 2A)]/(a + 3£). This formula, 
with A = 1.5 nm and £ = 9.7 nm as the best fit, actu- 
ally improves the agreement with the data for a < 50 
nm, where the deviation from the linear dependence is 
observable (see Fig. [5|. In our system, therefore, a fluid 
slip might be relevant for a < 50 nm (with a relatively 
large slip lengh £ ~ 10 nm). 

Summary. — We have developed a micro-manupulation 
technique allowing us to invert and amplify the move- 
ment of colloid particles due to thermophoresis of poly- 
mers induced by laser focusing. The polymer concentra- 
tion gradient created by thermophoresis causes migra- 
tion of a particle with a speed determined by the balance 
between the driving force and the viscous force. This 
new method and its application to trapping of molecules 
based on nonequilibrium effects does not rely on a specific 
character of particles and polymers, and thus provides 
further applications for manipulating a diverse range of 
colloidal particles as well as biological cells and DNA 
molecules [2TJ [30]. Our technique is also applicable with 
other water-soluble molecules [31]. We also note that the 
effect is stronger for PEG with a larger molecular weight 
at the same monomer concentration as far as it is lower 
than the overlap concentration. 
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